# Full text of "Masses of tetraquarks with two heavy quarks in the relativistic quark model"

## See other formats

HU-EP-07/21 HEPHY-PUB 842/07 Masses of tetraquarks with two heavy quarks in the relativistic quark model D. Ebert Institut fiir Physik, Humboldt-Universitdt zu Berlin, Newtonstr. 15, D- 12489 Berlin, Germany R. N. Faustov and V. O. Galkin Institut fiir Physik, Humboldt-Universitdt zu Berlin, Newtonstr. 15, D- 12489 Berlin, Germany and Dorodnicyn Gomputing Gentre, Russian Academy of Sciences, Vavilov Str. 40, 119991 Moscow, Russia W. Lucha Institute for High Energy Physics, Austrian Academy of Sciences, Nikolsdorfergasse 18, A-1050, Vienna, Austria Masses of tetraquarks with two heavy quarks and open charm and bottom are cal- culated in the framework of the diquark-antidiquark picture in the relativistic quark model. All model parameters were regarded as fixed by previous considerations of various properties of mesons and baryons. The light quarks and diquarks are treated completely relativistically. The c quark is assumed to be heavy enough to make the diquark configurations dominating. The diquarks are considered not to be point- like but to have an internal structure which is taken into account by the calculated diquark form factor entering the diquark-gluon interaction. It is found that all the {cc){qq') tetraquarks have masses above the thresholds for decays into open charm mesons. Only the I(J^) = 0(1+) state of {bb){ud) lies below the BB* threshold and is predicted to be narrow. PACS numbers: 12.40.Yx, 14.40.Gx, 12.39.Ki I. INTRODUCTION Recent experimental studies of the heavy meson spectroscopy revealed several new states, such as X(3872), F(4260), D!q(2317) etc., which cannot be simply accommodated in the quark- ant iquark (qq) picture [l[ . These states can be considered as indications of the possible existence of exotic multiquark states which were proposed long ago, e.g. in The idea to revisit the multiquark picture using diquarks has been put forward by Jaffe and Wilczek [3|. At present, vast experimental and theoretical evidence of the important role played by diquark correlations in hadrons is collected 2 The simplest multiquark system is a tetraquark, composed of two quarks and two anti- quarks. Heavy tetraquarks are of particular interest, since the presence of a heavy quark increases the binding energy of the bound system and, as a result, the possibility that such tetraquarks will have masses below the thresholds for decays to mesons with open heavy flavour. If such strong decays are kinematically forbidden, then corresponding tetraquarks can decay only weakly or electromagnetically and thus they should have a small decay width. In this paper we consider tetraquarks with two heavy quarks as bound systems of a diquark and antidiquark. Therefore we assume that both b and c quarks are heavy enough to make the attractive QQ interaction stronger than the Qq one. For tetraquarks contain- ing c quarks the obtained results crucially depend on this assumption. In particular, the doubly heavy {QQ'){qq') tetraquark {Q = b,c and q = u,d,s) is considered as the bound system of the heavy diquark {QQ') and light antidiquark {qq'), while the {cq){bq') tetraquark is the bound state of the heavy- light diquark (eg) and antidiquark {bq'). Masses of heavy tetraquarks with hidden charm {cq){cq) and bottom {bq){bq) were calculated in our previ- ous paper [5]. There the dynamical analysis has shown that X(3872) and y(4260) can be indeed the diquark-antidiquark tetraquarks with hidden charm. It was also argued that the corresponding ground-state tetraquarks with hidden bottom have masses below the open bottom threshold, and thus they should be narrow states. It is important to investigate the possible stability of the {QQ'){qq') tetraquarks since they are explicitly exotic states with the heavy flavour number equal to 2. Thus, their observation would be a direct proof of the existence of the multiquark states. Estimates of the production rates of such tetraquarks indicate that they could be produced and detected at present (SELEX, Tevatron, RHIC) and future facilities (LHC, LHCb, ALICE) Q. To calculate the masses of heavy tetraquarks we use the relativistic quark model based on the quasipotential approach in quantum field theory. Previously we considered in our model the mass spectra of the ground-state and excited doubly- heavy {QQq) [2] and heavy {qqQ) jsl, 0] baryons in the heavy-diquark-light-quark and light-diquark-heavy-quark approximations, respectively. The light quarks and light diquarks were treated completely relativistically. The internal structure of the light and heavy diquarks was taken into account by calculating diquark-gluon form factors in terms of the obtained diquark wave functions. Such scheme proved to be very effective and successful in our calculation of the masses of heavy baryons in good agreement with experimental data [lo|. The predicted masses of the T^b, and Sfc baryons proved to be very close to the recently measured ones [HI, H Q. Moreover, in Ref. [Qj it was shown that currently available experimental data on excited charmed baryons can be accommodated in the picture treating a heavy baryon as the bound system of the light diquark and heavy quark, experiencing orbital and radial excitations. This gives us additional confidence in the reliability of the diquark approximation within our model and motivates the consideration of tetraquarks as diquark-antidiquark bound systems. It is important to note that all parameters of our model were fixed in the previous calculations of meson mass spectra and decays, and we will keep their values in the following analysis of heavy tetraquarks. II. RELATIVISTIC QUARK MODEL In the quasipotential approach and diquark-antidiquark picture of heavy tetraquarks the interaction of two quarks in a diquark and the diquark-antidiquark interaction in a tetraquark are described by the diquark wave function (\l/rf) of the bound quark-quark state 3 and by the tetraquark wave function of the bound diquark-antidiquark state, respec- tively, which satisfy the quasipotential equation of the Schrodinger type [3] where the relativistic reduced mass is _ E,E, _ - (ml - ml)^ and El, E2 are given by IVP -ml + ml M"^ -ml + ml E-\ = , Eo = . (3) Here, M = Ei + E2 the bound-state mass (diquark or tetraquark), mi 2 are the masses of quarks (gi and ^2) which form the diquark or of the diquark {d) and antiquark {d') which form the heavy tetraquark (T), and p is their relative momentum. In the center-of-mass system the relative momentum squared on mass shell reads 2 ^ [M^ - (mi + m^YWM^ - (mi - m^f] The kernel l^(p, q; M) in Eq. ([1]) is the quasipotential operator of the quark-quark or diquark-antidiquark interaction. It is constructed with the help of the off-mass-shell scatter- ing amplitude, projected onto the positive-energy states. In the following analysis we closely follow the similar construction of the quark-antiquark interaction in mesons which were ex- tensively studied in our relativistic quark model [IJ, |15| . For the quark-quark interaction in a diquark we use the relation Vqq = Vgg/2 arising under the assumption of an octet structure of the interaction from the difference in the qq and qq colour states.^ An important role in this construction is played by the Lorentz structure of the confining interaction. In our analysis of mesons, while constructing the quasipotential of the quark-antiquark interaction, we assumed that the effective interaction is the sum of the usual one-gluon exchange term and a mixture of long-range vector and scalar linear confining potentials, where the vector confining potential contains the Pauli terms. We use the same conventions for the con- struction of the quark-quark and diquark-antidiquark interactions in the tetraquark. The quasipotential is then defined as follows 0, 15 (a) For the quark-quark {qq'), (Qq), {QQ') interactions, V{p,q; M) reads V{p, q; M) = Ui{p)u2{~p)V{p, q; M)ui{q)u2{-q), (5) with V(p,q;M) = i ^asD,Ml^l2 + Krnf(k)r^(k)r2;,(-k) + Kif(k) Here, is the QCD coupling constant; D^^, is the gluon propagator in the Coulomb gauge, D''{k) = -^, D^^{k) = -^(6^^-^], D°^ = D^° = 0, (6) ^ Obviously, it is important to study diquark correlations in gauge-invariant color-singlet hadron states on the lattice. 4 and k = p — q; 7^ and u{p) are the Dirac matrices and spinors u\p) e{p) + m \ Hp) 1 cr ■ p X (7) e{p) + m with e{p) = Vp^ + m'^. The effective long-range vector vertex of the quark is defined 3] by r^(k) = 7^ + ^(^,.uk' k = {0,k), (8) where k is the Pauh interaction constant characterizing the anomalous chromomagnetic moment of quarks. In configuration space the vector and scalar confining potentials in the nonrelativistic limit reduce to (l-£)Konf(r), with Ar + B, confv / ' '^confl where e is the mixing coefficient. (b) For the diquark-antidiquark {dd!) interaction, l^(p, q; M) is given by ^d^d <J ^yi^d'^d' +rAPWd'{n \Jd,,r,,v^oA^) + Kif(k)l MQ)MQ') (9) (10) where {d{P)\J fj\d{Q)) is the vertex of the diquark-gluon interaction which takes into account the finite size of the diquark and is discussed below P*^'^ = (-E^o, ±p) and Q^'^ = {E^o, ±q), Ed = (M2 - M| + MJ)/(2M) and Ed' = (M^ - M| + M|)/(2M) " . The diquark state in the confining part of the diquark-antidiquark quasipotential ffTTl) is described by the wave functions '1 for a scalar diquark, ^ Edip) for an axial- vector diquark. (12) where the four-vector £d{p) = [Sd ■ P) Md ,Sd + [^d ■ P)P Md{Ed{p) + Md) ^d{p)Pi^ = 0, (13) is the polarization vector of the axial- vector diquark with momentum p, Ed{p) = y P^ + M^, and erf(O) = (0, e^) is the polarization vector in the diquark rest frame. The effective long- range vector vertex of the diquark can be presented in the form J, 2\ EdEd d;ti ^JEdEd 2Md for a scalar diquark, for an axial-vector diquark. (14) 5 TABLE I: Masses M and form factor parameters (for definitions see Eq. (jl9p ) of light diquarks. S and A denote scalar and axial-vector diquarks, antisymmetric [q, q'] and symmetric {q, q'} in flavour, respectively. Quark content Diquark type M (MeV) ^ (GeV) C (GeV^) [u,d] S 710 1.09 0.185 {n,4 A 909 1.185 0.365 [u,s] S 948 1.23 0.225 {m,s} a 1069 1.15 0.325 {s,s} A 1203 1.13 0.280 where k = (0,k). Here, the antisymmetric tensor T,'^ is defined by iJ:,^)l = -tig,,5:-g,j;), (15) and the axial-vector diquark spin is given by {Sd-k)ii = —iskiu l^d is the total chromo- magnetic moment of the axial-vector diquark. The constituent quark masses rrih = 4.88 GeV, rric = 1.55 GeV, = = 0.33 GeV, rus = 0.5 GeV and the parameters of the linear potential A = 0.18 GeV^ and B = —0.3 GeV have the values typical in quark models. The value of the mixing coefficient of vector and scalar confining potentials e = —1 has been determined from the consideration of charmonium radiative decays [IJ] and the heavy-quark expansion [16|. The universal Pauli interaction constant k = —1 has been fixed from the analysis of the fine splitting of heavy quarkonia ^Pj- states jl3]- In this case, the long-range chromomagnetic interaction of quarks vanishes in accordance with the flux-tube model. III. LIGHT, HEAVY-LIGHT AND HEAVY DIQUARKS As the first step, we calculate the masses and form factors of the diquarks. As is well known, the light quarks are highly relativistic, which makes the v/c expansion inapplicable and thus a completely relativistic treatment is required. To achieve this goal in describing light and heavy-light diquarks, we closely follow our recent consideration of the spectra of light mesons and adopt the same procedure to make the relativistic quark potential local by replacing ei_2(p) = \l^i^2 + by £'1^2 (see the discussion in Ref. The resulting quark-quark interaction potential is equal to 1/2 of the qq interaction in the meson. We solve numerically the quasipotential equation with this complete relativistic potential which depends on the diquark mass in a complicated, highly nonlinear way. The obtained ground- state masses of scalar and axial- vector light and heavy diquarks [HI, 0, 0] are presented in Tables [T] and [m In order to determine the diquark interaction with the gluon field, which takes into ac- count the diquark structure, it is necessary to calculate the corresponding matrix element of the quark current between diquark states. This (diagonal) matrix element can be param- eterized by the following set of elastic form factors: (a) scalar diquark (5) {S{P)\J,\S{Q)) = h+{k'){P + Q),, (16) 6 TABLE II: Masses M and form factor parameters (for the definitions see Eq. (jl9|) ) of heavy- light and heavy diquarks. S and A denote scalar and axial-vector diquarks, antisymmetric [Q, q] and symmetric {Q, q} in flavour, respectively. 1 1 1 1 O ILdL K content \ It /"i 1 1 o y \^ Uii^lLdbL is. type o — — c D — h M (MeV) e (GeV) C (GeV2) M (MeV) e (GeV) C (GeV2) [Q,u] S 1973 2.55 0.63 5359 6.10 0.55 {Q,u} A 2036 2.51 0.45 5381 6.05 0.35 [Q,s] S 2091 2.15 1.05 5462 5.70 0.35 {Q,s} A 2158 2.12 0.99 5482 5.65 0.27 [Q,c] S 6519 1.50 0.59 {Q,c} A 3226 1.30 0.42 6526 1.50 0.59 {Q,b} A 6526 1.50 0.59 9778 1.30 1.60 (b) axial-vector diquark (A) {A{P)\J,\A{Q)) = -[el{P).e,{Q)]h{e){P + Q), +h,{e) {[em ■ Q]e,,,{Q) + [e,{Q) ■ P]4,(P)} +h3{k')^Js:iP) ■ Q][e,{Q) ■ P]{P + g)„ (17) where k = P — Q, and Sd{P) is the polarization vector of the axial- vector diquark (fT3|) . Using the quasipotential approach with the impulse approximation for the vertex function of the quark-gluon interaction, we find [1] hik') = 0, Ea + Mj (27r r^d P + 2e2(p) eiip + k) + ei{p) P Ed + Md J \ ti{p + k)+mi [2^ei{p + k)ei{p) vl/rf(p) + (1 ^ 2), (18) F(k^) 2yei(p + k)ei{p){ei{p) + mi). where are the diquark wave functions. We calculated the corresponding form factors F{r)/r, which are the Fourier transforms of F(k^)/k^, using the diquark wave functions found by numerically solving the quasipotential equation. Our estimates show that this form factor can be approximated with high accuracy by the expression F(r) = 1 - e (19) which agrees with previously used approximations [7|. The values of the parameters ^ and C. for light, heavy-light and heavy scalar diquark [q, q'] and axial- vector diquark {q, q'} ground states are given in Tables HI and HTl 7 IV. MASSES OF HEAVY TETRAQUARKS As the second step, we calculate the masses of heavy tetraquarks considered as bound states of diquark and antidiquark. For the potential of the diquark-antidiquark interaction (fTTD we eet ^ V{r) = Vboui(r) + Konf(r) + -< M2{E2 + M2) VL conf \ + Ei{Ei + Mi) E2{E2 + M2 + 2 1 + 1 \ K'o^nf(r) Ml M| >L- (S1 + S2) ^2(^2 + M2 1 \ V!^ Ml L- (S Ml (El + Ml Mi(Ei + Mi) M2{E2 + M2) E1E2 P Vcouiir) + K conf \ 1 + - r Ud 4 +- Ei_ Ml 4 VMi E^- ^^'^ M2 (Si - S2 -v;^oui(^) V^c'oul(0 + ^EiE2_ 4 M1M2 »+^Coul( (S1 + S2) r) 2r confl'") ^onf( X + - :Si-r)(S2-r)-Si-S2 Si • S2 (20) where 4 Fi(r)F2(r) 3 r is the Coulomb-like one-gluon exchange potential which takes into account the finite sizes of the diquark and antidiquark through corresponding form factors Fi 2(r). Here, Si 2 and L are the spin operators of diquark and antidiquark and the operator of the relative orbital angular momentum. Since we limit our considerations to the ground states of heavy diquark- antidiquark bound systems ((L^) = 0), the spin-orbit and tensor terms in the potential (120|) do not contribute in the further analysis. In the following we choose the total chromomag- netic moment of the axial- vector diquark = 0. Such a choice appears to be natural, since ^ In our paper [5| first two spin-orbit terms were missed. However they do not influence published numerical results, since masses mostly of ground states were calculated. Orbital excitations were considered only for the tetraquarks composed of the scalar diquark and scalar antidiquark for which the missed terms vanish. 8 TABLE III: Masses M of heavy-diquark (QQO^light-antidiquark {qq) states. T is the lowest threshold for decays into two heavy-light (Qq) mesons and A = M—T. All values are given in MeV. System State Q = Q' = c Q = Q' = b Q = c, Q' = b M T A M T A M T A iQQ'){ud) 0(0+) 7239 7144 95 0(1+) 3935 3871 64 10502 10604 -102 7246 7190 56 1(1+) 7403 7190 213 1(0+) 4056 3729 327 10648 10558 90 7383 7144 239 1(1+) 4079 3871 208 10657 10604 53 7396 7190 206 1(2+) 4118 4014 104 10673 10650 23 7422 7332 90 iQQ'){us) 7444 7232 212 4143 3975 168 10706 10693 13 7451 7277 174 7555 7277 278 4221 3833 388 10802 10649 153 7540 7232 308 4239 3975 264 10809 10693 116 7552 7277 275 ^(2+) 4271 4119 152 10823 10742 81 7572 7420 152 iQQ')iss) 0(1+) 7684 7381 303 0(0+) 4359 3936 423 10932 10739 193 7673 7336 337 0(1+) 4375 4080 295 10939 10786 153 7683 7381 302 0(2+) 4402 4224 178 10950 10833 117 7701 7525 176 the long-range chromomagnetic interaction of diquarks proportional to jjd then also vanishes in accordance with the flux-tube model. We substitute the quasipotential fl2U]) in the quasipotential equation ([T]) and solve the resulting differential equation numerically. The calculated masses M of tetraquarks with open charm and/or bottom composed from the heavy diquark, containing two heavy quarks {QQ', Q = b,c), and the light antidiquark {qq', q = u,d,s) are presented in Table UTTl In this table we give the values of the lowest thresholds T for decays into two corresponding heavy-light mesons [{Qq) = D^*\ Di*\ B^*\Bi*^], which were calculated using the measured masses of these mesons [lo[. We also show values of the difference of the tetraquark and threshold masses A = M — T. If this quantity is negative, then the tetraquark lies below the threshold of the decay into mesons with open flavour and thus should be a narrow state which can be detected experimentally. The states with small positive values of A could be also observed as resonances, since their decay rates will be suppressed by the phase space. All other states are expected to be very broad and thus unobservable. We find that the only tetraquark which lies considerably below threshold is the 0(1+) state of {bb){ud). All other {QQ'){qq') tetraquarks are predicted to lie either close to or significantly above corresponding thresholds. It is evident from the results presented in Table IIIII that the heavy tetraquarks have increasing chances to be below the open flavour threshold and thus have a narrow width with the increase of the ratio of the heavy diquark mass to the light antidiquark mass. 9 TABLE IV: Masses M of diquark (cg')-antidiquark [hq) states. T is the lowest threshold for decays into two heavy-light {Qq) mesons and A = M — T; T' is the threshold for decays into the Sc*^ and a light meson [q'q) and A' = M - T'. All values are given in MeV. System State q' = u q' = s M T A r A' M T A r A' icq'Mhu) 0+ 7177 7144 33 6818 359 7294 7232 62 6768 526 1+ 7198 7190 8 6880 318 7317 7277 40 6820 497 1+ 7242 7190 52 6880 362 7362 7277 85 6820 542 0+ 7221 7144 77 6818 403 7343 7232 111 6768 575 1+ 7242 7190 52 6880 362 7364 7277 87 6820 544 2+ 7288 7332 -44 7125 163 7406 7420 -14 7228 178 {cq'){bs) 0+ 7282 7247 35 6768 514 7398 7336 62 6818 580 1+ 7302 7293 9 6820 482 7418 7381 37 6880 538 1+ 7346 7293 53 6820 526 7465 7381 84 6880 585 0+ 7325 7247 78 6768 557 7445 7336 109 6818 627 1+ 7345 7293 52 6820 525 7465 7381 84 6880 585 2+ 7389 7437 -48 7228 161 7506 7525 -19 7352 154 In Table HVl the calculated masses M of tetraquarks composed from a (eg) diquark and a {bq) antidiquark are listed. ^ We also give the lowest thresholds T for decays into heavy-light mesons as well as thresholds T' for decays into the S^*-* and light (q'q) mesons and A'- '' = ]\^ _ )_4 gj-^^ ^j^g^i- Qj^iy 2+ states of {cq')(bq) have negative values of A and thus they should be stable with respect to decays into heavy-light {B and D) mesons. The predicted masses of lowest l"*" states of {cu)(bu) and {cu)(bs) tetraquarks lie only slightly above the corresponding thresholds T. However, all {cq){bq) tetraquarks are found to be significantly above the thresholds T' for decays into the -B^*-* and light (q'q) mesons. Nevertheless, the wave function of the spatially extended (cqMbq) tetraquark would have little overlap with the wave function of the compact Be meson 1^, thus substantially suppressing the decay rate in this channel. Therefore the above-mentioned {cq){bq) tetraquark states which are below the BD threshold have good chances to be rather narrow and could be detected experimentally. It is important to note that the comparison of the masses of heavy tetraquarks given in Tables UTTl and ITVl with our previous predictions [5| shows that the {QQ'){qq') states are, in general, heavier than the corresponding {Qq){Q'q') ones. This result has the following explanation. Although the relation [13] Mqq + < 2MQq holds between diquark masses, the binding energy in the heavy-light diquark (Qq)- heavy-light antidiquark (Qq) bound ^ Such (cm) (bu) tetraquarks were recently argued [18| to be the best candidates for experimental detection. ^ For the non-strange {cq)(bq) tetraquarks we give thresholds T' for decays of the / = states into B^*^ and T] or lu. These states should be more stable than the I — 1 ones, since their decays to Be and tt violate isospin. 10 TABLE V: Comparison of different theoretical predictions for the masses of heavy {QQ'){qq') tetraquarks (in MeV). System I{JP) this work [21] [22] [23] m [25] [26] [27] (cc) (ud) 0(1+) 3935 3931 3876 3764 3927 3905 4000 ± 200 1(0+) 4056 4150 1(1+) 4079 4186 1(2+) 4118 4211 {bb){ud) 0(1+) 10502 10525 10504 10558 10261 10426 10200 ± 300 1(0+) 10648 10587 10766 10690 1(1+) 10657 10712 10644 10774 10698 1(2+) 10673 10735 10790 10707 {bb)ius) 10706 10823 10680 10816 {bc){ud) 0(0+) 0(1+) 1(2+) 7239 7246 7422 7206 7244 7422 {be) {us) ^(2-^) 7572 7496 system is significantly larger than in the corresponding heayy diquark (QQ)-light antidi- quark {qq) one. This fact is well known from the meson spectroscopy, where heavy quarkonia QQ are more tightly bound than heavy-light mesons Qq. For instance, we found that some of the {cu){cu) tetraquarks lie below open charm thresholds while all ground-state {cc){ud) tetraquarks are found to be above such thresholds. In Table |V] we confront our results for masses of heavy {QQ'){qq') tetraquarks with other theoretical predictions [2l|, [22, l23|, M, l25|, l26|, l27| . In Ref. [21|] the authors solve the four-body problem using the expansion in the harmonic-oscillator basis in the framework of the nonrelativistic quark model with a phenomenological potential. The same model with a different expansion basis, which can accommodate asymptotic states of two free mesons, is applied for the calculation of heavy-tetraquark properties in Ref. [13]. The stability of tetraquarks with heavy flavours is studied by using a variational approach and tetraquarks are analyzed a nonrelativistic potential model in Ref. j23|. In Refs. [2^ in the chiral constituent quark model using a variational approach. The potential of this model includes one-gluon, confinement and meson-exchange interactions. The existence of a virtual tetraquark state ccud is discussed in Ref. [26[ on the basis of semi-empirical mass relations. QCD sum rules are applied for the QQud tetraquarks in Ref. [23]. The main difference between our approach and the above quoted papers consists in that, from the very beginning, we explicitly reduce the relativistic four-body problem to the subsequent solution of two relativistic two-body problems assuming the diquark-antidiquark structure 11 of the {QQ'){qq') tetraquarks. From Table IVl we see that most of the presented approaches predict that only the l"*" state of the (bb) {ud) tetraquark lies below the open-bottom threshold (see also (Isl). In Ref. [22[ it is claimed that also the l"*" state of the {cc){ud) tetraquark is weakly bound against the DD* threshold, if it has a molecular structure. Note that such structures are absent in our approach. A large binding energy in the 1^ state of the {cc){ud) tetraquark is found only in Ref. [13] and it is claimed to originate from the meson-exchange part of the quark interaction potential. The recent QCD sum rule analysis [23] finds that only the {bb)(ud) tetraquark is expected to be a narrow state. V. CONCLUSIONS In this paper we have calculated the masses of the ground states of tetraquarks with two heavy quarks assuming the diquark-antidiquark structure. Such approximation allowed us to reduce the very complicated relativistic four-body problem to the solution of two — significantly more simple — relativistic two-body problems. All considerations were done in the framework of the relativistic quark model which proved to be successful in describing numerous properties of mesons and baryons [l, SB [o, 14, 15, 16, 17]. The parameters of the model were fixed previously from the meson sector and are kept unchanged in the present analysis. The diquarks were treated dynamically. Their masses and form factors, which take into account the diquark structure, were calculated on the basis of a numerical solution of the corresponding relativistic quasipotential equation. Note that they are the same as in our previous studies of light and heavy diquarks in heavy [sl and doubly-heavy baryons 0], respectively. Light quarks and diquarks were treated completely relativistically without applying the v/c expansion. It was found that both the relativistic dynamics of light diquarks as well as their internal structure play an important role in the description of diquark- antidiquark bound states. The binding of a heavy diquark and a light antidiquark turned out to be weaker than the binding of a corresponding heavy-light diquark and heavy-light antidiquark. Thus, in contrast to the (eg) (eg) tetraquarks, which were discussed previously jsl, all the (ee)(gg') tetraquarks are predicted to be above the decay threshold into the open charm mesons. Only the I{J^) = 0(1+) state of {bb){ud) was found to lie below the BB* threshold. As a result, this state can decay only weakly and thus it should be narrow. The strange partner of this state {bb){us) is predicted to lie slightly (13 MeV) above the B*Bs threshold and, in principle, could be observed as a not too broad resonance decaying mainly into B*Bs. The investigation of the decay widths of heavy tetraquarks which are predicted to lie below the threshold of the open flavours represents another very important and interesting problem. It could be considered by means of the relativistic generalization of the analysis performed in Ref. 2^ . However, this problem is beyond the scope of the present paper and will be considered elsewhere. The authors are grateful to V. A. Matveev, D. Melikhov, M. Miiller-Preussker and V. I. Savrin for support and useful discussions. Two of us (R.N.F. and V.O.G.) were supported in part by the Deutsche Forschungsgemeinschaftmadei contract Eb 139/2-4 and by the Russian 12 Foundation for Basic Research under Grant No. 05-02-16243. For a recent review, see E. S. Swanson, Phys. Rep. 429, 243 (2006), and references therein. R. L. Jaffe, Phys. Rev. D 15, 267 (1977); Phys. Rev. Lett. 38, 195 (1977); V. A. Matveev and P. Sorba, Lett. Nuovo Cim. 20, 443 (1977). R. L. JafFe and F. Wilczek, Phys. Rev. Lett. 91, 232003 (2003). For a recent review, see R. L. Jaffe, Phys. Rep. 409, 1 (2005), and references therein. D. Ebert, R. N. Faustov and V. O. Galkin, Phys. Lett. B 634, 214 (2006). A. Del Fabbro, D. Jane, M. Rosina and D. Treleani, Phys. Rev. D 71, 014008 (2005). D. Ebert, R. N. Faustov, V. O. Galkin and A. P. Martynenko, Phys. Rev. D 66, 014008 (2002). D. Ebert, R. N. Faustov and V. O. Galkin, Phys. Rev. D 72, 034026 (2005). D. Ebert, R. N. Faustov and V. O. Galkin, arXiv:0705.2957 [hep-ph]. Particle Data Group, W.-M. Yao et al, J. Phys. G 33, 1 (2006). B. Aubert et al. [BaBar Cohaboration], Phys. Rev. Lett. 97, 232001 (2006). I. V. Gorelov [CDF Cohaboration] , |arXiv:hep-ex/ 0701056 . V. M. Abazov et al. [DO Collaboration], Phys. Rev. Lett. 99, 052001 (2007); T. Aaltonen et al. [CDF Collaboration], Phys. Rev. Lett. 99, 052002 (2007). D. Ebert, R. N. Faustov and V. O. Galkin, Phys. Rev. D 67, 014027 (2003). D. Ebert, V. O. Galkin and R. N. Faustov, Phys. Rev. D 57, 5663 (1998); 59, 019902(E) (1999). R. N. Faustov and V. O. Galkin, Z. Phys. C 66, 119 (1995); D. Ebert, R. N. Faustov and V. O. Galkin, Phys. Rev. D 73, 094002 (2006). D. Ebert, R. N. Faustov and V. O. Galkin, Eur. Phys. J. C 47, 745 (2006); Mod. Phys. Lett. A 20, 1887 (2005). H. J. Lipkin, |arXiv:hep- ph/0703190 , N. A. Tornqvist, Nuovo Cim. A 107, 2471 (1994). S. Nussinov and M. A. Lampert, Phys. Rep. 362, 193 (2002). B. Silvestre-Brac and C. Semay, Z. Phys. C 57, 273 (1993). D. Jane and M. Rosina, Few-Body Syst. 35, 175 (2004). D. M. Brink and Fl. Stancu, Phys. Rev. D 57, 6778 (1998). J. Vijande, F. Fernandez, A. Valcarce and B. Silvestre-Brac, Eur. Phys. J. A 19, 383 (2004). J. Vijande, A. Valcarce and K. Tsushima, Phys. Rev. D 74, 054018 (2006). B. A. Gelman and S. Nussinov, Phys. Lett. B 551, 296 (2003). F. S. Navarra, M. Nielsen and S. H. Lee, Phys. Lett. B 649, 166 (2007). H. J. Lipkin, Phys. Lett. B 172, 242 (1986); Phys. Lett. B 580, 50 (2004); M. Karliner and H. J. Lipkin, Phys. Lett. B 638, 221 (2006). [29] D. Melikhov and B. Stech, Phys. Rev. D 74, 034022 (2006).